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ABSTRACT 



Context. It has been proposed that the envelopes of luminous stars may be subject to substantial radius inflation. The peculiar structure 
of such inflated envelopes, with an almost void, radiatively dominated region beneath a thin, dense shell could mean that many in 
reality compact stars are hidden below inflated envelopes, displaying much lower effective temperatures. The inflation effect has been 
discussed in relation to the radius problem of Wolf-Rayet (WR) stars, but has yet failed to explain the large observed radii of Galactic 
WR stars. 

Aims. We wish to obtain a physical perspective of the inflation effect, and study the consequences for the radii of Wolf-Rayet (WR) 
stars, and luminous blue variables (LBVs). For WR stars the observed radii are up to an order of magnitude larger than predicted by 
theory, whilst S Doradus-type LBVs are subject to humongous radius variations, which remain as yet ill-explained. 
Methods. We use a dual approach to investigate the envelope inflation, based on numerical models for stars near the Eddington limit, 
and a new analytic formalism to describe the effect. An additional new aspect is that we take the effect of density inhomogeneities 
(clumping) within the outer stellar envelopes into account. 

Results. Due to the effect of clumping we are able to bring the observed WR radii in agreement with theory. Based on our new 
formalism, we find that the radial inflation is a function of a dimensionless parameter W, which largely depends on the topology of 
the Fe-opacity peak, i.e., on material properties. For W > 1, we discover an instability limit, for which the stellar envelope becomes 
gravitationally unbound, i.e. there no longer exists a static solution. Within this framework we are also able to explain the S Doradus- 
type instabilities for LBVs like AG Car, with a possible triggering due to changes in stellar rotation. 

Conclusions. The stellar effective temperatures in the upper Hertzsprung-Russell (HR) diagram are potentially strongly affected by 
the inflation effect. This may have particularly strong effects on the evolved massive LBV and WR stars just prior to their final 
collapse, as the progenitors of supernovae (SNe) Ibc, SNe II, and long-duration gamma-ray bursts (long GRBs). 

Key words, stars: Wolf-Rayet - stars: early-type - stars: variables: S Doradus - stars: interiors 



1. Introduction 

In the standard picture for the evolution of the most massive 
stars (with M ^ 30 M ), O stars have been proposed to evolve 
through successive luminous blue variable (LBV) and Wolf- 
Rayet (WR) phases befo re exploding as hydrogen (H) free su- 
pernovae (SNe) Ibc (e.g ., lLanger et aDll994t iMevnet & Maederi 
l2003l:lHeger et alj|2003l) . This evolutionary path is however still 
open as LBVs have more recently been sugge sted to be the 
direct progenitors of some H-rich type II SNe (|Kotak & Vinkl 
l2006tlSmith eTaT]l2007l iGal- Yam & LeonarJ2009l) . Additional 
renewed interest in WR stars arises from thei r connection to 
long-duration gamma ray burs ts (long GRBs, cf. lYoon & Langerl 
l2005tlWooslev & Hegerll2006h . 

What LBVs and WR stars have in common is their prox- 
imity to the Eddington limit. Not only i s this thought to be in- 
strumental for their mass-lo s s behavior (IVink & de Kotei 2002t 
Grafener & Hamarrnl l2008h IVink et all |201 It iGrafener et alj 
201 ll) . but it may also have key consequen ces for their stel - 



lar structures. Using mod ern OPAL opacities, Ishii et all dl999t) . 
iPetrovicet ail d2006l) , and lYungelson et al.1 d2008l) studied the in- 
terior structure of massive stars, revealing a "core-halo" config- 
uration that involves a relatively small convective core and an 
extended radiative envelope. This is referred to as the "inflation" 
of the outer envelope. 



It has been known for many decades that the observed radii 
of WR stars are almost an order of magnitude larger than canon- 
ical WR stellar structure models indicate. This is oftentimes at- 
tributed to their so-called pseudo-photospheres, which concern 
an "effective" photosphere several times larger than t he hydro- 
static radius, as a result of a dense stellar outflow (e.g. JCrowfherl 
l2007h . Pseudo-photospheres have also been discussed in the 
context of LBVs, although the issue is still under debate (e.g., 
ISmith et alJl2004 . 

An alternative explanation for the large radii of WR stars 
may involve envelope inflation. The issue is particularly relevant 
in the context of WR stars as GRB progenitors, as there have 
been several suggestions that the required radii of GRB pro- 
genitors are up to an order of magnitude ( ~ 1 vs. ~ 10^o) 
smaller than the radii of observed WR stars dModiaz et alj|2009t 
ICui et al.ll2010l) . The WR radius issue is thus highly relevant for 
GRB modeling. 

For LBVs the issue of their radii is equally interesting. The 
defining property of LBVs concerns their "S Doradus" cycles. 
On timescales of years to decades LBVs are seen to vary be- 
tween effective temperatures of ~30 kK (early B spectral type) 
to ~8 kK (early F). Due to the lack of convincing counter- 
evidence, these S Dor excursions in the upper Hertzsprung- 
Russell diagram have generally been assumed to occur at con- 
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stant bolometric luminosity (iHumphrevs & Davidsonfl994l) . but 
recen t studies have challenged this behavio r for a couple of ob- 
jects ( lGrohetal.H2009bt IClarket al.1 [20091) . Either way, the is- 
sue of the increased LBV radii during S Dor cycles remains 
un-questi oned, but a satisfying explanation for it has yet to be 
provided jVinkll2009l) . Outer envelope inflation may turn out to 
be an interesting explanation for S Dor variations, as will be de- 
tailed in the following. 

The paper is organized as follows. In Sect. [2] we present nu- 
merical stellar structure models that show the envelope inflation 
effect. In Sect.[3]we describe the effect analytically, including a 
new instability limit. We also provide a recipe that can be used 
to estimate the inflation effect for arbitrary stars with given (ob- 
served) parameters. The implications of our results are discussed 
in Sect.|U with focus on WR stars and LB Vs. The conclusions 
are summarized in Sect.[5] 

2. Model computations 

In this section we present stellar structure models for chemically 
homogeneous stars close to the Eddington limit. To model the 
conditions in classical WN-stars (i.e., WR stars in the phase of 
core He-burning that have lost their H-rich envelope), we use 
pure He-models, while models with hydrogen resemble the con- 
ditions in massive WNh stars, and LBVs. At this point we note 
that we only use the assumption of chemical homogeneity, to 
keep the models as simple as possible. The occurrence of an en- 
velope inflation is not connected to this assumption, as it only 
affects the outermost layers of the star. 

In Sect. l2.ll we give an overview of our numerical approach 
for computing the stellar structure, and in Sect. l2.2l we describe 
the envelope inflation, as it occurs in our models. In Sect. 12. 31 we 
present a model grid to investigate the occurrence of this effect 
in the HR diagram. Finally, in Sect. l2.4l we discuss the important 
influence of density inhomogeneities on our models. 

2. 1 . Numerical method 

Our models are computed with a code that integrates the stel- 
lar structure equations using a shooting m ethod. The numeri- 
cal me thods are described in the textbook by lHansen & Kawalerl 
(fl99l . The code is based on an example program that is dis- 
tributed with the boo k, but is completely re-wr itten, and updated 
with OPAL opacities dlglesias & Rogers1ll996l) . As the chemical 
structure of the stars is fixed, a set of four differential equations 
needs to be solved. These are the equations of hydrostatic equi- 
librium 



dP_ 

dm 



Gm 
"4^4' 



mass conservation 
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dr 

dm 4nr 2 p ' 

energy conservation 

dh _ 
dm 

and energy transport 
d\nT 



d\nP 



= V, 



(1) 



(2) 



(3) 



(4) 



pure stellar structure models, i.e. they include no time-dependent 
terms, the energy generation rate e in Eq.|3]does not take changes 
in the gravitational energy, due to contraction/expansion of the 
star into account. This omission has no consequences for the 
physics of the low-density envelopes discussed in this work, but 
can potentially affect the L/M ratio of WR stars in late evolution- 
ary stages (see Sect. l3.5.2l for a detailed discussion). 

Furthermore, we neglect dynamic terms due to mass loss, 
i.e., we only concentrate on the hydrostatic cas e. The dynamic 
case has been discussed bv lPetrovic et al. ( 20061) . For the case of 
WR stars with strong winds, iPetrovic et al.l found that dynamic 
effects may inhibit the formation of inflated envelopes. We will 
address this point in more detail in Sect. 13.5.31 

The choice of the outer boundary condition, particularly the 
temperature at the outer boundary, turns out to be important for 
the formation of inflated envelopes. Here, we prescribe p, and T 
at the outer boundary of the star in a way, that resembles reason- 
able values for the sonic points of stars with strong, radiatively 
driven winds. 

An envelope inflation only occurs in our models if we 
choose outer boundary temperatures below ~ 70 kK, well be- 
low the temperature of the hot Fe-opacity peak. Such temper- 
atures are expected at the sonic points of optically thick, radia- 
tively driven winds, such as the late-type WN stars described by 
IGrafener & Hamannl (120081) . The density p is chosen such that 
the resulting mass loss rate, 



M = 4npaRl, 



(5) 



where V = V ra d for the case of purely radiative energy transport, 
and V la d > V > V a d for the convective case. As our models are 



lies in a realistic range. 

We note that our envelope solutions are very robust with re- 
spect to the detailed choice of p, and T, as long as the tempera- 
ture T lies clearly below the regime of the hot Fe-opacity peak, 
which has its maximum at ~ 160kK. The effect of changes 
in the outer boundary condition is discussed in more detail in 
Sect.EO 

The stellar temperatures given in this work, are effective 
temperatures related to the outer boundary radius of our mod- 
els, i.e., T+ = (L/4-7TO-R+) 1 / 4 . For static atmospheres these val- 
ues will be almost identical to classical effective temperatures 
T e ff. In the presence of dense, extended stellar winds, our 
values will resemble the "effective core temperatures" T*, that 
are typically inferred fro m models for expanding atmospheres 
(cf. IGrafener et al.ll2002l) . 

2.2. Stellar structure models with inflated envelopes 

As a typical example for a model with an inflated envelope, we 
present a homogeneous stellar structure model for a 23 M Q he- 
lium star with solar metallicity (X = Y = 0.98 , Z - 0.02). 
In agreement with p revious works by llshii et ail dl999t) . and 
IPetrovic et al.l (120061) . this model develops an extended low- 
density envelope, with a density inversion close to the surface 
(cf. Fig.Q]). The star thus forms a shell around the actual stellar 
core, with a radius that is 3-4 times larger than the core radius. 

Fig. [2] shows that the inflated region is dominated by radi- 
ation pressure. In this model, the gas pressure P gas only ac- 
counts for a fraction of ~4-10~ 4 of the total pressure P. In con- 
trast to this, f gas and P ra d are roughly equal in the stellar core. 
The core is thus only fairly close to the Eddington limit, as re- 
flected by an Eddington parameter of T e « 0.4 (cf. Table[TJ. 
Here T e denotes the "classical" Eddington parameter, which is 
related to the Thomson opacity K e due to free electrons, i.e., 
F e = K e L rd d/(4-ncGM). The F e given in Tables[T]and|2] are com- 
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Fig. 5. Envelope solution in the P m d-P g as plane. The colors indi cate the logarithm of the Eddington factor F, for given P rac j and P gas , 
according to the OPAL opacity tables (Igl esias & Ro gers 1996), for our pure He model with 23 M . The envelope solution almost 
precisely follows a path with F = 1 . The radii of 2, 4, and 6 R e within the inflated envelope are indicated. Arrows indicate the slope 
of the solution, as expected from Eq.[T2l 




2 4 6 

r/R 

Fig. 1. Density vs. radius for our 23 M He model (cf. TableQ]). 



puted for a fully ionized plasma with a given hydrogen mass 
fraction X, and L la d = L. T e then simplifies to 

log(r e ) = -4.813 + log(l + X) + log(L/L ) - log(M/M ). (6) 

In contrast, the total Eddington factor T = KL m il{AncGM), 
which is related to the total opacity k, approaches unity in the 
inflated zone. This is shown in Fig. [3] where we compare the 
total luminosity within our model with the Eddington luminos- 
ity LecM = 4-ttcGM/k (the condition L ra d = LEdd is thus equiva- 
lent to T = 1). In Fig. [3] we can identify three zones within the 
star. In the inner, convective core, the total luminosity L is larger 
than the Eddington luminosity Z^dd- The main reason for this is 
that the energy transport is almost entirely convective, and thus 
£rad is very low. Above this region follows the radiative enve- 
lope with L < LEdd- On top of this follows the inflated zone, 
where LEdd almost equals L. T is thus extremely close to one, 
almost throughout the complete outer envelope of the star. This 
is surprising, because T is a function of the opacity k, which is a 
strongly varying function of p, and T. The conditions in the in- 
flated envelope thus demand for a mechanism that adjusts T very 
precisely to one. 
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2 4 6 

r/R Q 

Fig. 2. Pressure vs. radius. Gas pressure P gas , radiation pressure 
P m d, and total pressure P tot = P gas + f ra< j. 
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Fig. 3. Total stellar luminosity L(r), compared to the Eddington 
luminosity Lead = 4-7tcGM/k. 

To achieve a large radial extension it is necessary that the 
density scale height H is comparable to the stellar core radius R c . 
For our example model we have R c /H — GM(1 - T)l(a 2 R c ) ss 
3-10 3 (l - F). From Fig.|2]we can infer that (1 - T) a P gas /P ~ 
4T0~ 4 (cf. Eq.[T5l). Due to the low ratio of gas pressure to radi- 
ation pressure, F is thus just close enough to one, to achieve a 
large radial extension. 

A qualitative idea about the reason for the envelope inflation 
can be obtained from Fig.|4j where we plot the gas pressure P gas , 
and the total pressure P, vs. the radiation pressure P a <x, This di- 
agram is particularly useful, as for stars close to the Eddington 
limit, P is of the same order of magnitude as P ra d, i.e., in a log- 
arithmic diagram the stellar structure just follows a s traight line 
(in ac cordance with the Eddington standard model ( lEddingtonl 
Il918l) for which P r3 d/P = const, is adopted). For P gas the sit- 
uation is more complicated. In the stellar core, P gas follows a 
similar relation as P. In the extended envelope, P gds however 
drops significantly, and rises again after a certain minimum pres- 
sure P m j n is reached. As demonstrated in Fig. [5] the region where 
f gas drops corresponds precisely to the location of the Fe-opacity 
peak in the P ra d-Pgas plane. 
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Fig. 4. Stellar structure in the Prad-^gas plane. Gas pressure P gas , 
radiation pressure P rac j, and total pressure P tot = P gas + P rac j are 
plotted vs. P ra d throughout the whole star. 



More explicitly, starting from the stellar core, the solution 
follows a straight line with P gas /-P ra d = const. This naturally 
leads into a region in the P la d-^ gas plane, where the opacity k in- 
creases significantly, due to the presence of the Fe-opacity peak. 
When k reaches K^di, this leads to T — > 1. To avoid a super- 
Eddington situation, with T > 1, either L, a d, or k need to be 
reduced. L la d could in principle be reduced by convection, how- 
ever, in our case the involved densities turn out to be too low, i.e., 
convection is inefficient (cf. Sect. l3.5.4b . The only possibility to 
reduce a:, is thus to lower the density. For this reason the density 
has to drop in our model, until the tip of the Fe-peak is reached 
(cf. Fig.[5j. At the same time, however, this leads to a situation 
where P gas <k P ra d, and thus to F — > 1 (because, as already men- 
tioned above (1 — T) ~ Psas/P)- This mechanism thus forces the 
solution to follow a path with r = 1, i.e., with K(p, T) = ke^. 

As P ra d, and P gas are functions of p, and T, we can express k 
in the form K(P gas , P m &). This is actually done in Fig. [5] where the 
colors indicate K(P gas , A-adV^Edd on a logarithmic scale, as ob - 
tained from the OPAL opacity tables (llglesias & Rogers1ll996l) . 
In this plot we can see that the solution within the low-density 
zone is completely dominated by the topology of the opacity 
k, i.e., by material properties. As soon as the solution reaches 
low densities close to the Fe-opacity peak, it just follows a 
path with F = 1, until a minimum density p m i„ is reached. 
This minimum is located at the tip of the Fe-opacity peak, typi- 
cally at P ra d ~ 10 6 dyncnT 2 , corresponding to a temperature of 
~ 150kK. After the opacity peak is passed, the solution still fol- 
lows r = 1, but now with increasing P gas , leading to a density 
inversion. 

To maintain the described envelope inflation, it is thus nec- 
essary that 1) the star is close enough to the Eddington limit, so 
that r — > 1, and 2) convection is inefficient. If the envelope in- 
flation occurs, its properties are determined by the topology of 
the Fe-opacity peak in the Prad-Pgas plane, i.e., by material prop- 
erties. In Sect. [3] this will help us to find an analytical descrip- 
tion of the inflation effect. It also implies a dependenc e on mate- 
rial pr operti es like the metallicity Z, as discussed by llshii et all 
(1 19991) . and iPetrovic et al.l (120061) . In Sect.l2i4l we will discuss 
an additional important effect, namely the influence of density 
inhomogeneities within the inflated zone, which is based on a 
similar principle. 



4 



G. Grafener et al.: Stellar envelope inflation near the Eddington limit 



Table 1. Model parameters for pure helium models. 
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13.0 


5.371 


0.273 


5.09 


4.98 


1.1 




1.8 
















10.24 








12.0 


5.307 


0.255 


5.08 


5.01 


1.0 




1.5 
















6.97 








10.0 


5.155 


0.216 


5.07 


5.03 


0.9 


_ 


1.1 


_ 


_ 


_ 




_ 




_ 


3.07 


_ 


_ 


_ 


8.0 


4.959 


0.172 


5.05 


5.03 


0.8 




0.9 
















1.33 








X = 0.00, D = 


16 


































14.0 


5.429 


0.289 


5.08 


4.65 


1.2 


5.9 


8.5 


3.97 


0.80 


1.010" 


7 


7.410" 


n 


1.5 


14.59 


0.80 


3.97 


4.7- 10 s 


13.5 


5.401 


0.281 


5.08 


4.80 


1.2 


3.5 


4.3 


1.97 


0.66 


8.110" 


8 


8.9- lO- 


i i 


1.4 


10.99 


0.66 


1.97 


4.0- 10 5 


13.0 


5.371 


0.273 


4.88 


4.87 


2.8 




3.0 
















0.56 








12.0 


5.307 


0.255 


5.08 


4.95 


1.0 




1.9 
















10.49 








10.0 


5.155 


0.216 


5.07 


5.01 


0.9 




1.2 
















4.75 








8.0 


4.959 


0.172 


5.05 


5.02 


0.8 




0.9 
















2.06 









Notes. Given are: the stellar mass M*, the stellar luminosity L*, the classical Eddington parameter T e (Eq.[6]l, the effective temperature T c related 
to the core radius R c , the effective temperature T t related to the stellar surface radius R+, the inner core radius R c (at the bottom of the inflated 
envelope), the outer envelope radius R s , the stellar surface radius R+, the envelope extension AR in units of R c , and R c , the mass AM of the outer 
shell (with R s < r < R+), the "minimum density" p rain at the point of the lowest gas pressure within the inflated envelope, the ratio f p = p/p rain 
(with p from Eqs. ll8|[T9l >. the ratio Q according to Eq.[27] the parameters W c & W e according to Eqs.[l9][2T] and the total pressure P c at the outer 
envelope radius R e . The values given in this table are inferred numerically from our computational models, and may be compared to the analytical 
results from Sect.[3] 



2.3. A grid of homogeneous stellar structure models with 
envelope inflation 

To assess the effect of the envelope inflation in the HR diagram, 
we have computed a grid of homogeneous stellar structure mod- 
els with hydrogen mass fractions X = 0.0, 0.2, 0.4, and 0.7, that 
covers a luminosity range between log(L/L ) = 5.0, and 6.5. 
The results are compiled in TablesQ] and [2] In Fig. [6] we show 
the position of our models in the HR diagram, together with the- 
oretical zero-age main sequences (ZAMS) that do not take the 
inflation effect into account. 

At low luminosities, our He-burni ng models (with X = 0.0) 
closely resemble the He-ZAMS from IL anger! ( Il989bl). and our 
H-burning models (X = 0.2-0.7) the ZAMS from lSchaller et al.l 



Jl992h . Towards higher luminosities, our models develop in- 
flated envelopes, i.e., they show lower T+. With increasing lu- 
minosity, T+ decreases rapidly until, above a certain luminosity, 
we do not find static solutions anymore. The reason for the oc- 
currence of this important limit will be explained in Sect. 13.21 

In our present models, the envelope inflation sets in for 
Eddington factors T e between 0.3 (X = 0.0), and 0.4 (X = 0.7), 
i.e., within a relatively narrow range. The fact that it sets in ear- 
lier for H-deficient stars is a consequence of the fact that the 
L/M ratios of these objects are higher (although K e is lower for 
lower H abundances). In Sects. I3~2l & !3.3l we will show that this 
behavior, as well as the described upper luminosity limit, are a 
consequence of the strength, and shape of the Fe-opacity peak in 
the P ra d-Pg as plane. 
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Table 2. Model parameters for chemically-homogeneous stars containing hydrogen. 



m* u r e 


iog(r c ) 


log(T*) 


R c 




R+ 


AR 


AR 


AM 


Pmin 


fp 


Q 


w c 






[M e ] [L e ] 


[K] 


[K] 


[Re] 


[Re] 


[Re] 


[Rc] 


[«e] 


[Me] 


[±] 










r dyn , 


X = 0.70, D = 1 



155.0 


6.418 


0.435 


4.70 


4.20 


21.3 


163.3 


217.6 


6.68 


0.87 


7.1 


10" 




1.0 


10" 




2.3 


14.27 


0.87 


6.68 


8.6-10 5 


150.0 


6.397 


0.427 


4.70 


4.38 


20.7 


80.0 


90.6 


2.86 


0.74 


3.2 


10" 


4 


1.3 


lo- 


1 


2.2 


11.65 


0.74 


2.86 


8.6- 10 5 


145.0 


6.375 


0.420 


4.71 


4.47 


19.9 


53.6 


58.0 


1.69 


0.63 


6.2 


lo- 


-5 


1.5 


ur 


10 


2.2 


10.08 


0.63 


1.69 


8.6- 10 5 


140.0 


6.351 


0.412 


4.71 


4.53 


19.2 


41.1 


43.6 


1.14 


0.53 


2.2 


ur 


-5 


1.9 


lo- 


10 


2.2 


8.62 


0.53 


1.14 


8.5- 10 s 


130.0 


6.301 


0.396 


4.71 


4.60 


17.9 


29.2 


30.2 


0.63 


0.39 


5.5 


10" 


-6 


2.9 


ur 


10 


2.3 


6.45 


0.39 


0.63 


8.0- 10 5 


120.0 


6.247 


0.378 


4.71 


4.63 


16.7 




23.9 




















5.46 








110.0 


6.186 


0.359 


4.71 


4.66 


15.7 




20.2 




















3.77 








100.0 


6.118 


0.337 


4.71 


4.67 


14.6 




17.6 




















2.58 








90.0 


6.041 


0.314 


4.70 


4.67 


13.7 




15.7 




















1.76 








80.0 


5.952 


0.288 


4.70 


4.67 


12.7 




14.1 




















1.20 








70.0 


5.848 


0.259 


4.69 


4.67 


11.7 




12.7 




















0.81 








60.0 


5.724 


0.227 


4.68 


4.66 


10.7 




11.5 




















0.55 








50.0 


5.569 


0.191 


4.66 


4.65 


9.6 




10.2 




















0.36 








40.0 


5.368 


0.150 


4.63 


4.63 


8.9 




8.9 




















0.10 








30.0 


5.087 


0.105 


4.59 


4.59 


7.6 




7.6 




















0.09 








20.0 


4.647 


0.057 


4.53 


4.53 


6.0 




6.1 




















0.09 









X = 0.40, D = 1 



80.0 


6.179 


0.400 


4.73 


4.29 


14.3 


86.8 


106.6 


5.07 


0.84 


8.9- 10~ 4 


1.610 


10 


2.6 


12.21 


0.84 


5.07 


8.1-10 5 


75.0 


6.135 


0.385 


4.73 


4.50 


13.4 


35.7 


38.4 


1.66 


0.62 


2.2- 10- 5 


2.3- 10- 


10 


2.7 


9.46 


0.62 


1.66 


7.7- 10 5 


70.0 


6.088 


0.370 


4.72 


4.59 


13.1 


23.8 


24.8 


0.82 


0.45 


4.3-10- 6 


3.410- 


10 


2.3 


5.43 


0.45 


0.82 


7.6- 10 5 


65.0 


6.036 


0.354 


4.73 


4.63 


12.3 




19.0 














4.73 








60.0 


5.979 


0.336 


4.73 


4.66 


11.5 




15.7 














3.38 








50.0 


5.844 


0.296 


4.72 


4.68 


10.1 




12.0 














1.69 








40.0 


5.671 


0.248 


4.71 


4.68 


8.8 




9.8 














0.83 








30.0 


5.431 


0.190 


4.68 


4.67 


7.4 




8.0 














0.39 








20.0 


5.054 


0.120 


4.63 


4.63 


6.2 




6.2 














0.09 








15.0 


4.756 


0.080 


4.59 


4.59 


5.3 




5.3 














0.09 









X = 0.20, D = 1 



46.0 


5.975 


0.372 


4.74 


4.37 


10.7 


50.9 


59.5 


3.75 


0.79 


1.610 4 


2.3 


10" 


10 


2.4 


7.97 


0.79 


3.75 


7.8 


10 5 


45.0 


5.960 


0.367 


4.74 


4.45 


10.7 


36.7 


40.8 


2.43 


0.71 


4.1-10" 5 


2.6 


10" 


10 


2.2 


6.56 


0.71 


2.43 


7.7 


10 5 


44.0 


5.944 


0.362 


4.74 


4.50 


10.4 


29.0 


31.5 


1.78 


0.64 


1.610 -5 


2.9 


10" 


10 


2.2 


6.00 


0.64 


1.78 


7.7 


10 5 


42.0 


5.911 


0.352 


4.74 


4.57 


10.0 


20.9 


22.1 


1.09 


0.52 


4.2- 10~ 6 


3.8 


10" 


10 


2.2 


4.86 


0.52 


1.09 


7.5 


10 5 


40.0 


5.877 


0.341 


4.74 


4.61 


9.6 


16.8 


17.4 


0.75 


0.43 


1.7-10- 6 


4.9 


10" 


10 


2.3 


4.10 


0.43 


0.75 


7.0 


10 5 


35.0 


5.780 


0.312 


4.74 


4.67 


8.6 




11.9 
















2.73 










30.0 


5.664 


0.279 


4.73 


4.69 


7.7 




9.4 
















1.53 










25.0 


5.520 


0.240 


4.72 


4.70 


6.8 




7.8 
















0.84 










20.0 


5.334 


0.195 


4.71 


4.69 


6.0 




6.5 
















0.45 










15.0 


5.073 


0.143 


4.68 


4.66 


5.1 




5.4 
















0.23 










10.0 


4.664 


0.084 


4.61 


4.61 


4.3 




4.3 
















0.09 











X = 0.36, D = 1 (AG Car) 



73.0 


6.149 


0.397 


4.72 


4.23 


14.3 


104.3 


136.8 


6.30 


0.86 


2.2- 10~ 3 


1.610 


10 


2.2 


9.68 


0.86 


6.30 


8.0- 10 s 


70.0 


6.120 


0.387 


4.72 


4.44 


13.8 


46.1 


51.3 


2.36 


0.70 


6.8-10~ 5 


2.0- 10" 


10 


2.2 


7.69 


0.70 


2.36 


8.0- 10 5 


65.0 


6.069 


0.371 


4.73 


4.56 


12.7 


25.6 


26.9 


1.02 


0.50 


6.3-10- 6 


3.110 


10 


2.3 


5.89 


0.50 


1.02 


7.7- 10 5 


60.0 


6.013 


0.353 


4.73 


4.62 


11.8 




19.2 














4.90 









Notes. See Table[T]for an explanation of the listed physical quantities. 



2.4. The influence of density inhomogeneities 

For our pure He models in TableQ] the described envelope ex- 
tension starts at luminosities of log(L/L Q ) « 5.5, and fully de- 
velops for log(L/L ) 5.8 (cf. Fig.[6]l. As we will discuss in 
Sect. 14.11 the largest part of the Galactic H-free Wolf-Rayet stars, 
however, seems to show inflated envelopes for luminosities just 
below this value. In the following we will show that the infla- 
tion in this parameter range can be explained by density inho- 
mogeneities. Such inhomogeneities are expected to arise in the 
oute r envelopes of WR stars, e.g. due to strange -mode instabili- 
ties (iGlatzel & Kaltschmidtll2002t lGlatzelll2008l) . 

If material is inhomogeneous, or clumped, the density within 
clumps is larger than the mean density. As a result, the opacity k 
has to be evaluated for an increased density, instead of the mean 
density (note that k denotes the mass absorption coefficient in 
cm 2 /g). Here we assume a constant clumping factor D within 



the inflated envelope, and evaluate x(p, T) for the enhanced den- 
sity D x p, instead of the mean density p. This is equivalent to 
assuming a medium with clumps of a constant density Dxp, and 
an inter-clump medium which is void. In this picture the clumps 
would thus only fill a fraction / = l/D of the total volume. 
The same definition of clumping factors D, and filling factors 
/ is comm only used for the modeling o f inhomogeneous stellar 
winds (e.g. lHamann & Koesterkdfl998h . 

In the clumping approximation used here, it is assumed that 
no radiative flux gets lost in between clumps, e.g. by shadow- 
ing/porosity effects. Such a situation could be achieved by opti- 
cally thin clumps, or by a clump geometry (such as large shells, 
or pancakes) that is not sensitive to shadowing effects. 

Following our argumentation from the previous section, 
large envelope extensions, i.e., large scale heights H, are reached 
for Eddington factors F that are very close to one. We have 
shown that this is achieved for the low values of P gas at the tip 
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o 
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5.5 



5.0 



120M ® 



40M <f 



. He-ZAMS 



20 M * 




X=0.7 



4.8 4.6 

log(T,/K) 



Fig. 6. HR diagram for pure helium models (red), and models 
containing hydrogen (blue). 



of the Fe-opacity peak, as P g!is /P ~ (1 - T). By the assump- 
tion of clumping, the opacity peak is shifted towards even lower 
mean densities, i.e., towards lower (mean) values of P gas . In this 
way, the envelope inflation is further enhanced (a more concise 
explanation of this effect is given in Sect. [3] Eqs. [T9ll20ll29T l. 

In Fig.[6]we show He-ZAMS models that are computed with 
clumping factors D — 4, and D = 16. As expected, the enve- 
lope extension occurs much earlier in these models. In Sect. 14.11 
we will show that these models cover the observed HR diagram 
(HRD) positions of the Galactic H-free WR stars, i.e., we can 
explain the observed temperatures of these objects with moder- 
ate clumping factors. The adopted values for D are in notable 
agreement with spectroscopically determined clumping factor s 
for the winds of WR stars, e.g. bv lHamann & Koesterkd(ll998l) . 



3. Analytical description of the envelope inflation 

In the previous section we have shown that the properties of in- 
flated envelopes are chiefly determined by the fact that the so- 
lution has to follow a path with F — 1 in the P m d-P ' gas plane 
(cf. Fig . O . In the following we elaborate on this, and derive an 
analytical description of this process. As we concentrate on the 
physics of the inflated envelopes alone, the resulting relations 
do not depend on the internal structure of the star, i.e., they are 
generally applicable to stars with given (observed) parameters. 
In Sect. 13.11 we start with the equations describing the envelope 
structure. In Sect. 13.21 we derive analytical expressions for the 
radial extension of the envelope AR, and in Sect. l3.3l we present 
a recipe to estimate AR based on observed/adopted stellar pa- 
rameters. Furthermore, we derive an estimate for the mass of 
the surrounding shell, AM, in Sect. 13.41 Finally, we discuss in 
Sect. 13.51 to which extent the inflation effect may be affected by 
model assumptions. 



3.1. Inflated low-density envelopes near the Eddington limit 

In this section we concentrate on the physics of the outermost 
stellar envelope, where m = M, and L = const. With these 
assumptions, we eliminate two equations from the stellar struc- 
ture equations (Eqs.QJO, so that only the equation of hydrostatic 



equilibrium (Eq.[TJ, and the energy transport equation (Eq.|4ji are 
left. 

In Sect. l2.2l we concluded, that two conditions are manda- 
tory to maintain an inflated envelope. 1) the star needs to be 
close enough to the Eddington limit, so that F = 1 can be 
reached, and 2) convective energy transport needs to be ineffi- 
cient. Under condition 2) (the validity of this condition is dis- 
cussed in Sect. l3.5.4l ). the energy transport equation (Eq.[4| be- 
comes purely radiative, i.e., L = L m &, and V = V ra d, with 



V rad = 



Pk L 



\6nac T 4 GM' 



(1) 



After substituting Eq.Q] and Eq.|2]in Eq.|U we thus obtain the 
energy transport equation in the form 



4.2 — = - 



dT 
57 



KpL 



4aT 3 4nr 2 c' 
If this relation is written in terms of P ra( j 
8P md 4-aT 3 8T k P L 



(8) 



(a/3)T we obtain 



dr 



3 dr 



Anr 2 c 



(9) 



for the (outward directed) radiative acceleration g md . Together 
with the equation of hydrostatic equilibrium 



~Pg> 



dP _ d(P gas + P rad ) _ GMp 
dr dr r 2 

we obtain by division of Eq.fTOland Eq.|9] 



dP dP„ 



dP mA dP rdd 



+ 1 



,'rad 



1 

r' 



or 



dPo 



dP 



rad 



i 

f 



L. 



(10) 



(ID 



(12) 



Notably, dP sas /dP m d depends only on the Eddington factor 
T, i.e., via r/F e = k/k s on T e , and k. This is in line with our pre- 
vious finding that the envelope inflation is fully determined by 
the classical Eddington factor F e (as given by Eq.|6]l, which is a 
function of stellar parameters, and k, which depends on material 
properties. 

In Fig. [5] we compare our numerical solution from Sect. 12. 21 
(blue line) with the slopes dP rll d/dP gd>i (black arrows) follow- 
ing from Eq.[T2] in the P ra d-Pgas plane. The background colors 
indicate the value of F = r e (/c//c e ) for our model, where k has 
been extracted from the OPAL opacity tables dlglesias & Rogers! 
fl996h . 

According to Eq.[12] the point at which the envelope solu- 
tion crosses the Eddington limit, i.e. where T = 1, needs to be 
an extremum in P gas . Fig.|5]nicely shows, that our numerical so- 
lution indeed approaches the Eddington limit, and crosses it at 
the point of the lowest possible gas pressure. If the Eddington 
limit would be crossed earlier, at a higher density, the gas pres- 
sure had to increase so fast, that extremely high densities would 
be required at the outer boundary. The only way to reach reason- 
ably low temperatures, and densities at the outer boundary, thus 
leads around the Fe-opacity peak, via a low gas pressure, i.e., 
via low densities. The formation of a low-density envelope with 
a density inversion is thus a natural consequence of the topology 
of the Fe opacity peak in the Aad-Pgas plane, in combination 
with Eq.frU 

For the example in Fig. [5] we find that P gas /P m d ~ 4-10~ 4 at 



the tip of the Fe-opacity peak. In Sect. l2.2l we already discussed 
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that this implies that F — > 1, as P gas /P « (1 - T). This rela- 
tion follows from Eq.QT] if we assume a slowly changing ratio 
Pmd/P- In this case we have <5(lnP rad )/<9(lnP) ~ 1 (cf. Fig.©. It 
thus follows that 



<9(lnP rad ) P 8P 



rad 



3(lnP) Aad OP 
which leads to 

or 
P 



1, 



ad 



P 



+ Pm 



(i-r). 



(13) 



(14) 



(15) 



As discussed in Sect. 12. 21 this small value of (1 - F) leads to a 
density scale height of the order of the stellar radius, i.e., to an 
envelope inflation. 

The occurrence of an envelope inflation is thus intimately 
connected to the fact that P gas <k P lad at the tip of the Fe-opacity 
peak. As noted earlier, the location of this point in the P ra d-Pgas 
plane is only determined by the topology of the opacity peak, 
i.e., by material properties. In the following we will take ad- 
vantage of this fact, to obtain analytical estimates of the radial 
envelope extension AR. 



3.2. Computation of the radial extension AR 

To compute the radial extension AR, we integrate Eq.|9] and 
Eq. [TO] from the outer boundary of the inflated zone, at the "en- 
velope radius" P e , to the inner boundary, at the "core radius" P c . 
As this integration is performed within the inflated zone, where 
P = Pi-ad, and k = keh, Eqs.|9]and[T0]degenerate to one equation, 
that describes the change of P rad , and thus also the change of the 
temperature, in hydrostatic equilibrium. 



<?P r ad K Edd pL GMp 



dr 



4nr 2 c 



(16) 



or 



(17) 



dPrad 

= GMd-. 

P r 

The integration is performed from R e , where P = P e , to P c , 
where P = P c , leading to 



GM 



r 



dp 



rad 



i 



Rc 



ip c p 

The left hand side of this equation can be written in the form 
R c AP 



(18) 



W c s 



GM p 



(19) 



where AP = P c — P e , and 1/p is the mean of 1 jp with respect to 
P ra d- In Sect. l3.3l we will show that AP, and p can be estimated 
alone on the basis of opacity tables. For given W c , the radial 
envelope extension can then be computed from Eq.Q~8] which 
becomes 



Re 

R, 



1 



1 - W c 



= 1 + We, 



(20) 



where the latter equality uses a definition of W in terms of the 
observable, outer-envelope radius R e 



W, = 



Pe AP 

GM p ' 



(21) 



The expression for the radial extension AP 
W c , or W e follows directly from Eq.[20l 



1 



AR _ 

P7 ~ l-W c 



I, 



P e - P c for given 



(22) 



AP 

pT 



= 1 - 



1 



1 + W e 



(23) 



For the W c , and W e as determined from our models (TablesQ] 
and[2]i these relations are fulfilled very precisely. 

The parameters W c / e are of essential importance for the in- 
flation effect. As u = 3P/p in a radiation dominated gas (where 
u denotes the specific energy per volume), AP/p is related to the 
specific energy per gram of material within the inflated zone. 
According to Eq.[18] it denotes the increase of the specific en- 
ergy, from the outer to the inner edge of the inflated envelope, 
divided by 3. The parameters W c / e indicate the ratio between 
this energy increase, and the gravitational energy per gram of 
material, at P c / e . 

Notably, Eq.[20]imposes an instability limit for W c — > 1, for 
which P e — > oo. This limit coincides with the limit that we have 
described in Sect. 12. 31 where we found that no hydrostatic solu- 
tions exist, if certain luminosities are exceeded. Stable solutions 
only exist for W c < 1, which is equivalent to 



AP GM 



(24) 



As AP/p is largely fixed by material properties, Eq.|24] im- 
poses a limit on the gravitational energy at the core radius P c . If 
the gravitational energy becomes too small (i.e., P c becomes too 
large), the stellar envelope becomes energetically unbound, and 
there exists no hydrostatic envelope solution. This limit might be 
of profound importance for the occurrence of instabilities close 
to the Eddington limit, as observed, e.g. for LB Vs (cf. Sec. l4.2t . 

The occurrence of the envelope inflation is thus mainly due 
to the fact that the Fe-opacity peak forces the envelope to low 
densities p, while the temperature of the opacity peak is fixed, 
i.e., u = const. This leads to a very high specific energy u/p, i.e., 
the envelope becomes less gravitationally bound. 

3.3. A recipe to compute AR 

The remaining step for the computation of AP is to obtain esti- 
mates of the parameter W (Eqs.[19j |2T1 i, which depends on the 
ratio M/R (where P denotes either P c , or P e ), and the ratio AP/p. 
While M, and P are given stellar parameters, AP/p needs to be 
determined from opacity tables. 
Based on Eqs.[18j[T9]we have 



AP 
P 



r 



dP 



rad 



(25) 



In Fig.|7]we illustrate that, due to the topology of the Fe-opacity 
peak, the integral on the right hand side is roughly equal to 
AP/(2p min ), i.e., 



W : 



P AP 



P AP 



GM p GM fpPmin 



(26) 



with f p 2. Numerically obtained values for f p from our model 
computations confirm this value (cf. Table [TJ and [2j. W can thus 
be computed on the basis of AP = P c - P e , and p m j n . 

To extract these values from the opacity tables, we make use 
the fact that the envelope solution follows almost precisely a path 
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/(10°dyn cm 2 ) 



Fig. 7. Plot of 1/p vs. f rac j. The dashed horizontal lines indicate 
1/Pmin, and l/(10p m in). The intersections of the latter with the 
envelope solution represent boundary pressure P c , and P e , as in- 
dicated by vertical dashed lines. The grey shaded area illustrates 
that the integral j dP/p ss AP/(2p min ). 



with r = 1 , i.e., a path with k = /CEdd in the P ra d-Pgas plane. 
From this path we extract the minimum gas pressure P m [ n , at the 
tip of the Fe-opacity peak (cf. Fig . O ■ We use the density at this 
point as our reference density p m j n . For given P mm , we define the 
boundaries of the inflated envelope as the two adjacent points on 
the path with k - where P gas = fp ■ P m ; n (in this work we 
use fp = 10). The radiation pressure at these two points gives us 
P c , and P e . Our choice of fp = 10 has a moderate influence on 
the resulting values, will however not change the overall results. 

The procedure described above thus provides p lrun , and AP = 
P c - P e for given K Edd , or T e (note that r e (/f Edd //e e ) = 1), i.e., we 
can now compute W, and AR. To determine AP/p mm over a wide 
parameter range, we employ a root finding algorithm to ext ract 
AP, and p min from the OPAL tables dldesias & Rogersll 19961) . as 
described above. The value of AP/p m j n depends on the chemical 
composition, i.e., X and Z, and on F e . 

At this point we define a function Q(X, Z, r e ), that only de- 
pends on APj p m i n 



Q(x,z,r e ) 



RqAP 
GM G p mi 



(27) 



Adopting a Galactic metallicity of Z = 0.02, we determine a 
fitting relation for Q(X, F e ), based on the numerically determined 
values of AP/p mm . The resulting relation has the form 



log(0 = -2.135- 0.617 X + 8.453 F e 



(28) 



The fits are displayed in Fig. [8] where the crosses denote the nu- 
merically determined values of Q, for the range X = 0-0.7, and 
the blue lines denote the fitting relation Eq.|28] for X — 0, and 
X = 0.7. 

The values of Q determined this way, lie between Q ~ 1, 
and Q ~ 150 (cf. Fig.[8}, with T e = 0.3-0.6. For lower Q (i.e., 
lower F e ), the extension of the Fe-peak becomes so small that it 
does not cover a factor fp = 10 in f gas . Under these conditions 
an envelope inflation may still occur, however only of moderate 
strength. For Q larger than ~ 150, p m j n becomes so small that 
the range of the opacity tables is exceeded. In this case we still 
expect a very strong envelope inflation, but are just not able to 
infer Q from the tables. 



a 

5 1 




Fig. 8. Fitting relation for Q(X,Z = Z Q ,r e ). Black crosses indi- 
cate numerically determined values of Q, for X = - 0.7. Blue 
lines indicate our fitting relation Eq.|28] for X — (top), and 
X = 0.7 (bottom). 



For given Q, W can be determined 
additionally allow for clumping with a 
(cf. Sect. l2.4t . we have p = (f p /D)p m i n , i.e.. 
the form 



from Eq.|26] If we 
clumping factor D 
we can express W in 



W = 



RAP 



_ Q(x,z,r e ) (r/Rq) 

GMp min (f p /D) <JJD) (M/M Q ) ' 



(29) 



with f p x2. 

To compute the envelope extension AR for a star with given 
parameters M, L, X, and R (where R may denote R c or R e ), the 
following steps thus need to be performed. 

1) For given M, L, and X, T e can be computed from Eq.|6] 

2) For given X, and r e , Q can be estimated from Eq.l28l 

3) For given Q, W can be determined from Eq.[29] adopting 
reasonable values for f p , and D. Note that the stability limit 
Eq.[24]can be expressed in the form Q/{f p /D) < M/R (in 
solar units). 

4) Finally, for given W, the radial extension AR can be com- 
puted from Eq.|22] or Eq.l23l 

E.g., for our He models in Table[T](with X = 0), masses are 
typically of the order of 15 M , and core radii are ~ 1 /? sun , i.e., 
M/R ~ 15 in solar units. This means that an unstable situation, 
with W = 1, is reached for Q/(f p /D) ~ 15 (cf. Eq.|29]l. For f p = 
2, and D = 1, this corresponds to Q ~ 30. According to Eq.l28l 
this value is reached for F e ~ 0.4, in agreement with our models 
with D = 1 in Table|2] For a large clumping factor D = 16, 
a value of W = 1 is already reached for Q ~ 2, suggesting 
an instability at F e ~ 0.3, in agreement with our models with 
D = 16 in Table|2] 

In Fig. [9] we c ompare our ana lytical predictions, based on 
core r adii R c from lLangerl (Il989bl) . with our numerical models. 
Long-dashed lines indicate predicted HRD positions within our 
analytical formalism for clumping factors D = 1,4, and 16, and 
an adopted f p =2. The short-dashed line indicates the results for 
fp = 1-5, a nd D = 1. Note that the use of the He-ZAMS models 
by lLangerl introduces inconsistencies. Nevertheless, our formal- 
ism describes the envelope inflation very well, in particular its 
dependence on the clumping factor D. Also the uncertainty due 
to the adopted f p turns out to have only moderate influence on 
the results. 
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Fig. 9. HR diagram for pure He models (red), compared with our 
analytical relations. Long-dashed lines: analytical relation with 
f p — 2, short-dashed line: analytical relation with f p = 1.5. 

3.4. The mass of the surrounding shell AM 

The mass AM of the dense shell surrounding the inflated zone 
(cf. Fig.QTJi may be of interest, e.g. in relation to the energy bud- 
get, and the timescales of variations in AR. 

To compute AM we use the following assumptions. 1) The 
radial extension of the shell is small compared to the total stellar 
radius, i.e., r = R e . 2) The mass of the shell is small compared 
to the stellar mass, i.e., m = M. 3) The pressure P at the stellar 
surface is small compared to the pressure at the outer boundary 
of the inflated region P e . 

Under these assumptions, the equation of hydrostatic equi- 
librium (Eq.[JJ) can be integrated directly, and we obtain 



AM = P e 



4< 
GM 



(30) 



The mass of the shell thus simply follows from the fact that the 
pressure force P e 47r/? 2 at the rim of the low density region needs 
to be in equilibrium with the gravitational force GMAMjR\ on 
the shell. Again, P e depends on material properties. Using P e « 
0.4...0.9T0 6 dyn/cm 2 from our models (cf. TableQ] and we 
obtain 



AM/M e * 0.5. ..1.0-10 



-9 (RJRq) 4 

M/M ■ 



(31) 



AM thus mainly depends on the radius of the star. This is in 
agreement with our model computations (TableQ] and [2]) where 
we find AM as low as 10~ 9 M for compact WR stars, whereas 
our larger, LBV- type models reach AM up to 10~ 2 M Q . 

3.5. Dependence on model assumptions 

The occurrence of an envelope inflation, involving a density in- 
version, has been reported in many previous w orks. Examples 
are the Wolf-Rayet models by ISchaererl d!996l) . and models for 
very massive m ain sequence stars, as well as hydrogen-deficient 
carbon stars, by ISaio et al 1 (119981) . What all these stars have in 
common, is their proximity to the Eddington limit. Especially 
the latter work, which discusses the connection of inflated en- 
velopes to the strange-mode instability, demonstrates the uni- 
versality of the principles discussed here, for stars covering the 
range of 0.9-1 10 M . 



o 
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Fig. 10. The effect of p and T at the outer boundary, on the enve- 
lope solution. The solid blue line indicates our standard 23 M Q 
He-model in the P ra d-Peas plane. The dashed lines indicate so- 
lutions for which the boundary values for the density p, and 
the temperature T have been changed. The inflated zone be- 
tween 2 and 6 R Q (as indicated by the radius marks) is mainly 
affected by the choice of the temperature T at the outer bound- 
ary. Temperatures T > 80 kK, can lead to a cut-off of the inflated 
zone. 



On the other hand, the inflation effect has not been detected 
in many standard works on stellar evolution . E.g., i n the rotating 
evolutionary models by iMevnet & Maederl d2003l) . the 120 M 
track enters the H-free WR stage with a mass of 23.2 M Q , and 
an effective core temperature of = 123 kK. According to our 
results in Table[TJ we would expect a considerable envelope in- 
flation, with T+ ~ 60 kK, for the same object. In this section 
we discuss possible sources of such discrepancies, such as the 
influence of the imposed outer boundary condition CSect. l3.5TTb . 
dynamical terms on a secular timescale (Sect. 13. 572b . dynamical 
terms introduced by mass-loss (Sect. 13. 53b . and the treatment of 
convection and turbulence (Sect. 13. 574b . 

3.5.1 . The effect of the outer boundary condition 

In Fig. [10] we illustrate the influence of the adopted boundary 
values for p, and T, on the envelope solution. The blue line 
indicates our standard 23 M Q He-model, with boundary values 
T = 50kK, and p = 5T0~ 10 g cm 2 . In addition we show six test 
models (dashed lines), for which we have changed p to values 
of 5T0~ 8 g ctrT 2 and 5T0~ 12 g cnT 2 , for temperatures of T — 50, 
80, and 160kK. 

Firstly, the envelope solution turns out to be very robust with 
respect to changes of the outer boundary. For increasing f ra d, 
i.e., towards the interior of the star, all test models merge into the 
standard solution. However, boundary temperatures in the region 
of the Fe-opacity peak, i.e., temperatures that would lie within 
the inflated zone, lead to a cut-off of the inflated envelope. In 
Fig.QJjJ the inflated zone is indicated by radius marks (2-6 R G ). 
Models with outer boundary temperatures T > 80 kK display en- 
velopes where the inflated zone ends close to the respective outer 
boundary temperature, i.e. the radii are significantly reduced. 

Boundary conditions that take the back-warming effect of 
an optically t hick w ind into account, as the ones discussed 
by ISchaerer] (fl996h . may thus inhibit an envelope inflation 
iScl 



Schaererl indeed detected a very similar effect in his model 
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computations). In the same way, strong stellar winds with 
high temperatures at the w ind base, as the ones described by 
iGrafener & Hamanr] d2005l) for early WC subtypes, will prevent 
the formation of inflated envelopes. This is not surprising, as 
such winds are initiated by the hot Fe-opacity peak, i.e., the Fe- 
opacity peak lies within the wind. We thus expect a strong influ- 
ence of the detailed wind physics on the inflation effect. 

3.5.2. Dynamical terms on a secular timescale 

The contraction/expansion of a star during its evolution can lead 
to energy sources/sinks, in addition to the nuclear energy pro- 
duction rate e, which is included in our models (Eq.0. To esti- 
mate the importance of this effect, we compare the gravitational 
energy E g , which is gained, or lost on a timescale t, with the stel- 
lar luminosity L, and estimate in this way a limiting timescale 



Table 3. Estimated limiting mass-loss rates Mq according to 
Eq.[33] based on an adopted value of p m j n = 5TCT 11 g/cm 3 . 



A£ g _ AM(M - AM)G 
~L~ ~ LR ' 



(32) 



Here AM is the mass of the layers involved in the contrac- 
tion/expansion, and R is a typical radius. 

Because Wolf-Rayet stars are very compact, and evolve on 
relatively short timescales, such gravitational terms may play a 
role in their evolution. However, as the masses of the inflated en- 
velopes discussed here are very small, these terms play no role 
for the inflation effect. For our 23.7 M He-model, which repre- 
sents the most extreme case of an inflated envelope, we estimate 
r = 0.5 d, based on the parameters from TableQ] and R - R c . 
We thus expect no influence of secular changes on the inflation 
effect. In Sect. l4.2l we will discuss the case of LBVs, for which 
gravitational terms indeed seem to play a role on timescales of 
the observed S Dor-type variability. 

For the complete radiative envelope of a typical WR star 
(here we use the 15 M model from TableQ] with AM = 1.5 M , 
and R = 0.4 R ) we estimate r = 5000 yr. Particularly in the end 
phases of the WR evolution, gravitational terms may thus play an 
important role in determining the L/M ratio of WR stars, which 
is an important input for our formalism. 

3.5.3. Dynamical terms due to mass-loss 

Because of the low densities, mass loss may introduce con- 
siderable velocity fields within the inflated envelopes, via the 
equation of continuity M = Anpvr 2 . The effect of the result- 
ing dyn amical terms on the en velope structure has been investi- 
gated bv lPetrovic et alJ (12006). They found that strong mass loss 
can inhibit the formation of inflated envelopes, if the velocities 
within the inflated zone exceed the local escape speed, i.e., when 
v(dv/dr) =s v 2 /r has the same order of magnitude as GM/r 2 . This 
condition imposes an upper limit for the mass-loss rate 



MG 

Mo=4n PrrAsi R J m \j— , 

Am 



(33) 



where we have taken into account that the maximum velocity 
is reached close to R m = (R e + R c )/2, where p w p m \ n . An in- 
spection of the dynamical term that would result from the den- 
sity structure within our own models, suggests however that this 
limit should be lowered by ~ 0.4 dex. 

In Table[3] we have compiled conservative estimates of 
Mo, based on p^n = 5T0~ n g/cm 3 , for a set of WR, and 
LBV models. For the WR models with the lowest tempera- 
tures, which may be subject to an envelope inflation, we find 



log(L) 


X 




M 




R* 
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log(Mo/M yr ') = -4. 2. ..-3.7. Fo r comparison, s pectro scopi- 
cally determined mass-loss rates bv lHamann et al 1 d2006h reach 
values of up to log(M/M yr~') = -4. 3. ..-4.1, dependent on L. 
In view of the involved uncertainties, it is thus not clear whether 
the envelope inflation of WR stars is affected by dynamical ef- 
fects. 



3.5.4. Convection and turbulence 



As already noted in Sects. 12721 and l3.ll the inefficiency of con- 
vection is a mandatory condition for the occurrence of an en- 
velope inflation. As the convective efficiency may be affected 
by details of the numerical implementation, or effects that are 
not taken into account in the standard mixing-length approach, 
we give basic physical arguments, why the convective efficiency 
should be low for the cases discussed in this work. 

Generally, the strong increase in opacity near the iron 
bump will lead to the onset of convection through the usual 
Schwarzschild stability criterion. But for the models here such 
convection should be quite inefficient in the sense that it can 
carry only a small fraction of the local stellar flux F. To see 
this, note that a robust upper limit to the flux carried by con- 
vection is set by the free-streaming of internal energy U at the 
local sound speed v s , Ff s = u^U. In general, U = t/ gas + f/ ra( j, but 
in the luminous stars here, and particularly in the region around 
the iron bump, the radiative energy dominates^. The associated 
maximum fraction of stellar flux F that could be carried by con- 
vection is thus 



F 



1 F 



— (f) 00034 



(34) 



where T c is the effective temperature for the (pre-inflated) core 
radius R c . The last equality applies the characteristic iron bump 
temperature T ~ 1.5 x 10 5 K, with associated sound speed 
v s « 50km/s. This shows directly that for core effective tempera- 
tures appropriate for Wolf-Rayet stars, convection can carry only 
a small fraction of the total stellar flux, viz. 0.3% for T c = 10 5 K, 
and about 5% for T c = 50, 000 K. This implies that, even when 
convection sets in, the bulk of the stellar flux is still transported 
by radiative diffusion, thus justifying this basic simplifying as- 
sumption in our analysis of envelope inflation. 

We note that, due to the dominance of radiation pressure, 
and radiation energy in the inflated zone, the contributions from 



1 Such convective transport of radiative energy assumes radiative dif- 
fusion from the convective cell takes longer than the turnover time. If 
this breaks down, the convective flux upper limit defined here would be 
even lower. 
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internal gas pressure/energy are generally negligible with re- 
spect to Ff S . The same holds for turbulent pressure/energy, as 
^turb = (2/3)wturb = (l/3)pa 2 M 2 is of the same order of mag- 
nitude as gas pressure/energy (cf.. lMaederll2009l p. 105). Here 
"turb denotes the turbulent energy density, a the sonic speed, and 
M the Mach number, which is expected to be lower than one. 

For low convective efficiency, the structure of inflated en- 
velopes is solely determined by Eq.Q~6] which is (in first order 
approximation) equivalent to the condition that F — 1 . The lat- 
ter only follows from material properties, so that the addition of 
turbulent pressure has almost no effect on the resulting envelope 
structure. 



4. Implications 

Let us next discuss the implications of the predicted envelope 
inflation for stars close to the Eddington limit, i.e., for WR stars 
and LBVs. Both types of objects are extremely luminous, and 
are thought to reside close to the Eddington limit. In Sect. 14. II 
we focus on the still unexplained WR locations in the HR di- 
agram. In Sect. 14. 21 we discuss the possible connection of the 
radius inflation to the S Doradus radius changes characteristic of 
LBVs. 



4. 1 . The radius problem for Wolf-Rayet stars 

In Fig.QTI we show an HR diagr am for galactic WR s tars, as 
obtained from spectral analyses bv lHamann et al.l ( I2006I) . As al- 
ready highlighted by these authors there exists a radius problem 
for H-free WR stars, namely that these objects are not located 
on the theoretical He- main sequence. A ccording to a population 
synthesis performed bv lHamann et al.l almost all WR stars with 
a hydrogen surface abundance X < 0.05 are expected to dis- 
play stellar temperatures in excess of lOOkK. However, in 
the observed HR diagram they cover a broad temperature range 
from 40-1 40 kK, with the majority below lOOkK (red symbols 
in Fig.rTTt. The observed WR ra dii are thus larg er than predicted 
by stellar structure models (e . g . . lLan gerll 1 989bf ) . 

At this po int we note that the stellar temperatures T*, as 
determined by lHamann et al 1 (120061) . are not classical effective 
temperatures (related to t — 2/3), but effective temperatures, 
related to the (hydrostatic) surface radius i.e., 



(35) 



star. Moreover, the strong temperature sensitivity of He-burning 
guarantees an almost fixed temperature in their centers. The 
main unknown for the determination of the stellar structure is 
thus the opacity. In fact, with the advent of the OPAL opacity 
data (llglesias & Rogers! 1996b and their strong Fe-opacity peaks, 



{4-ncr SB R 2 



with the Stefan-Boltzmann constant cr SB . R+ is the inner bound- 
ary radius of the employed atmosphere models, and is typically 
located at a much higher optical depth of r max * 20. As long as 
R+ is located in the hydrostatic part of the wind, its value does 
no t depend on the detailed choice of r max Q. The determined 
bv lHamann et al.l thus reflect actual surface radii, which are not 
affected by the optical depth of the strong stellar winds of WR 
stars. The low in Fig.Qj] thus display a true inconsistency 
with respect to the stellar structure models. 

This is particularly surprising, as H-free WR stars are very 
simple objects. Due to their large convective cores, and high 
mass-loss r ates they are ex pected to be close to chemical ho- 
mogeneity dLangerlll989al) . Even if the material in their He- 
burning cores is partly processed to carbon and oxygen, the 
mean molecular weight stays almost constant throughout the 



llshii et al j d 1 999b were the first to predict the large envelope in- 



2 The cases where R+ is located above the hydrostatic surface, are 
indicated as upper limits in Fig.ll II 



flation that we also describe in this work. 

llshii et all (1 19991) could however not explain the HRD posi- 
tions of these stars. In the galaxy, most H-free WR stars are lo- 
cated at luminosities below l og(L/L Q ) » 5.8, reaching down to 
log(L/L ) » 5.3 (cf. Fig.rm>. l lshii et al] demonstrated that unre- 
alistically high metallicities, up to Z = 0.1, would be needed to 
achieve a significant radius inflation in this regime. An important 
result of the present work is that the same effect can be achieved 
by the assumption of an inhomogeneous (clumped) structure 
within the inflated zone (cf. Sect l2.4t . For moderate clumping 
factors, up to D = 16, our models cover the observed parame- 
ter range (red (dashed) lines in Figs.l6landfTTI). Notably, similar 
clumping factors are determined spectroscopic ally, for the winds 
of WR stars (e.g jHamann & Koesterkelll998l) . 

Th e m etallicity dependenc e, as predicted by llshii et all 
(1 19991) and lPetrovic et al.l (120061) . may however still play an im- 
portant role in the determination of the temperatures of WR 
stars, as their su btype distribution seems to depend on Z (e.g., 
ICrowtherl 120071) . Later spectral subtypes are typically found at 
higher Z, which is in line with a stronger radius inflation due to 
the stronger Fe-opacity peak at higher Z. This effect may how- 
ever also depend on the in creased mass loss at higher Z (see 
ICrowther & Hadfieldll2006l) . We furt her note that, as disc ussed 
in Sect. 13.5.31 the mass-loss limit by iPetrovic et alj (|20Q6:) may 
inhibit an envelope inflation for WR stars with strong mass-loss. 

4.2. Radius changes of S Doradus-type LBVs 

A defining property of LBVs are S Doradus-type radius vari- 
ations on timescales of months to years. The inflation effect 
provides a reasonable explanation for radius variations on such 
timescales, as only small portions of the stellar envelope are af- 
fected. The un derlying reason for th e time dependence may be 
manifold. E.g.. IGeorgiev et aTl d20 lib recently proposed that the 
variability of HD 5980, an LB V-like star that temporarily shows 
a WR-type spectrum, may be triggered by a transition from a 
hot, to a cool wind base. This picture fits very nicely into the 
framework of the inflation effect. As discussed in Sect. 13. 5. Tl 
we expect that a hot, optically thick wind would inhibit a radial 
inflation, i.e., a wind transition as proposed by IGeorgiev et al.l 
could indeed lead to substantial radius changes. 

In the following we discuss the case of AG Car, one of 
the be st-studied galacti c S Doradus LBV s (e.g . jLeithereret al.l 
[19941: IStahl et al.fedoTt IVink & de Koteil 120021) . More recently, 
AG Car has been st udied in detail throughout its full S Dor cy- 
cle (1985-2003), bv lGroh et al.l d2009bll20Tll) . During this pe- 
riod the star increased its radius from 58.5 R e to 120.4 R Q , cor- 
responding to a change in the stellar temperature from = 
26, 450 K to 16,650 K. Such radius and temperature changes 
are characteristic aspects of the S Dor cycle, but they remain 
as y et unexplained by stell ar structure calculations (see how- 
ever IStoAers^Chia [l993|). For AG Car the luminosity was 
found to decrease from log(L/L ) = 6.17 to 6.0 between 
the S Dor visual m inimum and maximum phase. Furthermore, 
iGroh et aTl (12009b) determined a considerable hydrogen defi- 
ciency for AG Car, with a hydrogen surface mass fraction of 
X = 0.36. 
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Fig. 11. Herzsprung-Russel diagram of the Gala ctic WN stars, and the LBV AG Car. Red/blue symbols indicate observed HRD 
positions of WN stars from lHamann et all (120061) (blue: with hydrogen (X > 0.05); red: hydrog en-free (X < 0.05)) . Black symbols 
indicate the HRD positions of AG Car throughout its S Dor Cycle from 1985-2003, according to lGroh et al.l (l2009bl) . Large symbols 
refer to stars with known distances from cluster/association membership. The symbol shapes indicate the spectral subtype (see 
inlet). Arrows indicate lower limits of for stars with strong mass loss. The observations are compared to stellar structure models 
with/without hydrogen from this work (blue/red lines, cf. Fig.|6]l, and for AG Car (X = 0.36, black line). 



We have computed a series of homogeneous stellar struc- 
ture models for the relevant parameter range, with a hydrogen 
mass fraction of X = 0.36, and masses ranging from 60-73 M G 
(cf. Table|2]i. These models display a similar luminosity range as 
AG Car, and partly show very low T*, down to ~ 17,000 K, due 
to a substantial radius inflation. A comparison with the observed 
HRD positions of AG Car is shown in Fig.Ql] Notably, during 
the minimum phases of the S Dor cycle, AG Car (indicated by 
black symbols) shows a very good agreement with our homoge- 
neous models (black dashed line), i.e., its position in the HR di- 
agram can be explained by a chemically homogeneous star with 
M ~ 70 M Q , and a substantial radius inflation. 



The core of such a star would be relatively small, with 
R c ~ 14 R Q , corresponding to an effective core temperature 
T c ~ 53 kK (cf. Table|2]i. Due to the envelope inflation, the star 
would display a core-halo structure, as described in Sect. 12.21 
According to Eq.jJTJ the mass AM, of the outer shell around the 
inflated envelope, is expected to change throughout the S Dor cy- 
cle. For the smallest radius (58.5 R G ), we obtain AM ~ 10~ 4 M Q , 
while AM ~ 2-lO~ 3 M for the largest extension (120.4 ^ Q ). The 
energy to lift this material from the core radius to the stellar sur- 
face is A£ g ~ AMMG/R C = 3.8-10 46 erg. This value compares 
well with the 'missing energy' due to the decrease i n luminos- 
ity AEj , lm = 1.83-10 39 erg/sx575d = 9.1-10 46 erg (cf. iGroh et all 
l2009bl) . We thus find that within the inflated envelope, the S Dor- 
type variability introduces dynamical effects as the ones dis- 
cussed in Sect. l3.5.2l 



The good agreement with our models suggests that the exten- 
sion of the envelope of AG Car may be explained by the inflation 
effect. The mass-l oss rate of AG Car (log(M/M yr _1 ) < -4.2, 
Groh et al. 2009b) lies well below the upper mass-loss limit dis- 
cussed in Sect l3.5.3l (cf. TableO, implying a stable, static con- 
figuration that is not affected by dynamical effects due to mass 
loss. Based on our chemically homogeneous models, AG Car 
lies close to the Eddington limit, with an Eddington factor of 
T e = 0.39. According to Eq.[28j this results in a Q parameter of 
~ 9, and with Eq.|29]we obtain W c ~ (9/f p )(RJR )/(M/M Q ) ~ 
0.8 (where we have adopted f p = 2.2 from Tabled . The star is 
thus very close to the instability limit (W c = 1), as discussed in 
Sect.ll2l 

One possible mechanism that could lead to an unstable situ- 
ation involving Wc > 1, is the reductio n of the effective stellar 
mass through rotation. iGrohet al.l (120 111) found that the observed 
rotational velocity of AG Car declines from u mt sin i = 220 km/s 
in the minimum phase, to 90 km/s for the largest extension. For 
these values the ratio of centrifugal to gravitational acceleration 
at the stellar surface decreases from v^ ol R+/(GM) = 0.21/ sin 2 ;, 
to 0.07/ sin 2 i. The corresponding reduction of the effective stel- 
lar mass by 20% in the minimum phase, could thus just lead to 
a situation where W c ~ 1 . In this case the outer stellar envelope 
would become unstable, and could start expanding. The reduc- 
tion of the rotational velocity during the expansion would lead 
to W c < 1, and a re-stabilization of the envelope. 

IGroh et all d2009al 1201 lb proposed a somewhat similar sce- 
nario for AG Car, where the Eddington limit is exceeded due 
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to rotation. In this work, we have demonstrated that already for 
moderate Eddington factors of T e ~ 0.43 AG Car may form an 
inflated stellar envelope that is dominated by radiation. Within 
this envelope the star adjusts very precisely to the Eddington 
limit (F = 1). In our scenario, the instability that potentially in- 
duces the S Dor cycle is due to the fact that the specific internal 
energy due to radiation exceeds the gravitational potential, with 
the stellar envelope thus becoming unbound. 

We want to point out that the formation of inflated envelopes 
does not depend on the assumption of chemical homogeneity in 
our numerical models. The instability follows from our formal- 
ism in Sect.[3j and is thus independent of the internal structure 
of the star. For AG Car, the observed luminosity (log(L/L Q ) = 
6.17), and surface hydrogen mass fraction (X = 0.36) imply 
a mass of 74 M fo r the case of chemical homogeneity (cf. 
iGrafener et al.ll201lh . In reality, AG Car is likely more chemi- 
cally enriched in the core, i.e., its mass is lower, and its core 
radius larger. This would lead to an even earlier onset of the in- 
flation effect, and its instability. As the evolutionary status of 
LBVs is not well known, it is however not possible to obtain 
reliable mass estimates for this type of object. 



5. Conclusions 

In the present work, we have discussed the possible formation 
of radially inflated stellar envelopes for sta rs approaching th e 
Edd ington limit, as orig inally proposed by llshii et all d!999t) . 
and lPetrovic et alj ( 12006b . In addition to numerical models, we 
could provide an analytic description of this process, leading 
to the discovery of a new instability limit, and a clumping de- 
pendence. Both effects turn out to be of profound importance. 
While the former may be connected to S Doradus-type instabil- 
ities in LBVs, the latter can account for the large observed radii 
of Galactic WR stars, and thus resolve the long standing WR 
radius problem. 

Within our analytical approach we could describe the infla- 
tion effect in terms of a dimensionless parameter W (Eqs. [28l[29b 
that characterizes the ratio between internal and gravitational en- 
ergy at the base of the envelope. W can be estimated on the basis 
of opacity tables and a combination of stellar parameters. The 
envelope inflation occurs when W approaches one. For W > 1, 
we find that the envelope becomes gravitationally unbound, i.e., 
there exists no static solution. A prerequisite for envelope in- 
flation is the proximity of the star to the Eddington limit, with 
Eddington factors of T e > 0.3-0.35, at solar metallicity. For a 
given chemical composition and F e , the resulting stellar radius 
depends on the ratio M/R, and on the degree of (in)homogeneity 
of the material within the envelope (characterized by a clumping 
factor D). 

Envelope inflation has a strong impact on the effective tem- 
peratures of stars in the upper HR diagram. For Wolf-Rayet stars 
it can account for the 'radius problem', i.e., that the observed 
stellar temperatures of H-free WR stars are much lower than 
predicted by canonical stellar structure models. To reproduce the 
observed HRD positions of these objects, it is necessary to as- 
sume that the material within the envelope is clumped. The re- 
sulting clumping factors lie in the range of D = 1-16, in notable 
agreement with t ypical clumping factors dete cted in the winds 
of WR stars (e.g jHamann & Koesterkdll998l) . As the envelope 
inflation is expected to depend on metallicity (llshii et alJll999t 



3 Note that, by definition, T c only takes the opacity due to free elec- 
trons into account. 



iPetrovic et al.l l2006h . it may also account for the observed Z- 
dependence of the spectral subtype distribution of WR stars, i.e., 
that WR stars i n high-Z environ ments show lower effective tem- 
peratures (e.g.. ICrowtherll2007h . 

For the luminous blue variable AG Car we suggest that the 
observed HRD position of a ~ 70 M star, with a relatively com- 
pact core, might be subject to substantial envelope inflation. 
During its S Dor cycle, AG Car is found to increase its radius by 
a factor of two, while its luminosity decreases by a factor of 1.5. 
The luminosity decrease might be explained as a result of the 
formation of a dense outer shell with a mass of - 2-lO~ 3 M as 
predicted by our models. The energy needed to lift this material 
from the core to the outer radius is of the same order of magni- 
tude as the 'missing energy' due to t he decrease in luminosity . 
Similar to the mechanism proposed bv lGroh et al.l j2009all201 lb . 
the S Dor variability of AG Car may be due to the effect of stellar 
rotation, which reduces the effective stellar mass, and thus leads 
to an instability with W > 1, that forces the star to expand. 

We conclude that the stellar effective temperatures in the up- 
per HR diagram are potentially strongly affected by the inflation 
effect. The peculiar structure of the inflated envelopes, with an 
almost void, radiatively dominated region, beneath a thin and 
dense shell could mean that many, rather very compact stars, are 
hidden below inflated envelopes, and thus display much lower 
effective temperatures to the observer. Also compact, fast rotat- 
ing stars may be hidden this way. This may particularly affect 
massive stars just before the final collapse, and the question of 
whether they form a SN Ibc, or a GRB (for the latt er smaller 
core r adii are inferred than the observed WR radii, cf. ICui et al] 
120101) . The complex envelope str ucture may also affect the early 
X-ray afterglow of GRBs (e.g. 120071) . 

Acknowledgements. We thank A. Maeder for a very constructive referee report, 
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